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Abstract 

rj meson production in both proton-proton and proton-neutron collisions is 
investigated within a relativistic meson exchange model of hadronic inter- 
actions. It is found that the available cross section data can be described 
equally well by either the vector or pseudoscalar meson exchange mechanism 
for exciting the 5ii(1535) resonance. It is shown that the analyzing power 
data can potentially be very useful in distinguishing these two scenarios for 
the excitaion of the 5ii(1535) resonance. 
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I. INTRODUCTION 



The production of rj mesons in nucleon-nucleon {NN) collisions near the threshold energy 
has been a subject of considerable interest in the past few years, since the existing data are by 
far the most accurate and complete among those for heavy meson production. Consequently, 
they provide an opportunity to investigate this reaction in much more detail than any of 
the other heavy meson production reactions. In addition to the total cross section for the 
pp — >• ppr] reaction we now have data for — pnr] 0] and — >• dr] [|,|| reactions. The 
differential cross section data for the pp ppr] reaction are also available. Consequently, 
there have been a large number of theoretical investigations on these reactions |]TU|-|I^. 

The production of rj mesons in NN collisions is thought to occur predominantly through 
the excitation (and de-excitation) of the 5*11(1535) resonance, to which the rj meson couples 
strongly. However, the excitation mechanism of this resonance is currently an open issue. 
For example, Batinic et al. have found that both vr and rj exchanges are the dominant 



excitation mechanisms. However, they have considered only the pp — * pp-q reaction. Faldt 
and Wilkin |]13| and Gedalin et al. fl^ have considered both the pp — > ppr] and pn — > pnr] 



reactions. In the analysis of Ref. |T3| the pn dr] reaction is also considered. These 
authors [p!3l , p!4| find p exchange to be the dominant excitation mechanism of the 5*11 (1535) 
resonance. In particular, it has been claimed |TB] that p meson exchange is important for 
explaining the observed shape of the angular distribution of the pp ppr] reaction. In an 
anaylsis of pp ppr] reaction, Santra and Jain also considered p meson exchange as the 
dominant excitation mechanism of the 5*11(1535) resonance. In contrast to the findings of 
Refs. [|T^-[l5[|, Peiia et al. have found that the dominant contribution arises, not from the 
5ii(1535) resonance current, but from the shorter range part of the nucleonic currents. In 
this work, we shall report on another possible scenario for exciting the S*ii(1535) resonance 
that reproduces both the pp — > ppr] and pn pnr] reactions and discuss the possibility of 
disentangling these reaction mechanisms. 

Although we focus here on the problem just mentioned, the description of r] meson 
production in NN collisions presents other interesting aspects. For example, the r] meson 
interacts much more strongly with the nucleon than do mesons like the pion so that not only 
the NN final state interaction (FSI), but also the r]N FSI is likely to play an important role, 
thereby offering an excellent opportunity to learn about the r]N interaction at low energies. 
In fact, the near-threshold energy dependence of the observed total cross section for r] meson 
production differs from that of vr and r]' production, which follow the energy dependence 
given simply by the available phase-space together with the NN FSI. The enhancement of 
the measured cross section at small excess energies in r] production compared to those in 
TT and r( production is generally attributed to the strong attractive r]N FSI. In addition to 
all of these issues, the theoretical understanding of r] meson production in NN collisions 
near threshold in free space is also required for investigating the dynamics of the S*ii(1535) 
resonance in the nuclear medium, the possible existence of r]NN bound states, and the 
possibility of using r] to reveal the properties of high-density nuclear matter created in 
relativistic heavy-ion collisions. 

In section II we introduce our model and define the meson production currents whose 
details are given in the appendix. An alternative model, which is similar to previous works 
based on p exchange dominance, is introduced in section HI. The results are given in section 
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IV. Section V provides a summary. 



II. THE MESON EXCHANGE MODEL 

Our model of the A^A^ —>■ NNrj reaction is based on a relativistic meson exchange model 
of hadronic interactions. The reaction amplitude is calculated in the Distorted Wave Born 
Approximation. Here we follow a (non-rigorous but otherwise economic) diagrammatic 
approach to present our formulation . A more rigorous derivation of the reaction amplitude 
will be reported elsewhere. We start by considering the meson-nucleon (MN) and NN 
interactions as the building blocks for constructing the total amplitude describing the NN — >■ 
NNM reaction. We then consider all possible combinations of these building blocks in a 
topologically distinct way, with two nucleons in the initial state and two nucleons plus a 
meson in the final state. In this process of constructing the total amplitude, care must be 
taken in order to avoid diagrams that lead to double counting. Specifically, the diagrams 
that lead to mass and vertex renormalizations must be excluded since we choose to use the 
physical masses and coupling constants. The resulting amplitude constructed in this way is 
displayed in Fig. |I[ The ellipsis indicates those diagrams that are more involved numerically 
(including, in particular, the MN FSI, which otherwise would be generated by solving the 
three-body Faddeev equation). So far there are very few attempts to account for them 
T^ , p!7[| . These higher order terms are also not considered in this work. 



In order to make use of the available potential models of NN scattering, we will carry 
out our calculation within a three-dimensional formulation which is deduced from Bethe- 
Salpeter formulation by restricting the propagating two nucleons to be on their mass shell. 
We follow the procedure of Blankenbecler and Sugar . The meson production amplitude 
illustrated in Fig.l then takes the following familar form 

M = (1 + Tj-^hG^f^*)J{l + 2G';+)t/+)) , (1) 

where T(^i,f) denotes the NN T-matrix interaction in the initial(i)/final(/) state, and G^ij) 
is the three-dimensional Blankenbecler-Sugar (BBS) propagator. The superscript ± in T(jj) 
as well as in in Eq.(|I]) indicates the boundary condition, (— ) for incoming and (+) 

for outgoing waves. The production current is denoted by J, which is defined by the MN 
T-matrix with one of the meson legs attached to a nucleon (first diagram on the r.h.s. in 

Fig. E), 

J = ^ [T{mn^m'n)] ^ iPhi' \^m'nn\2 + (1^2), (2) 

M' 

where T(_mn^m'n) stands for the MN T-matrix describing the transition M'N MN, 
^M'NN and Pm' stand for the M'NN vertex and the corresponding meson propagator re- 
spectively. The subscripts 1 and 2 stand for the two interacting nucleons 1 and 2. The 
summation runs over the intermediate meson M'. Eq.{^ is the basic formula on which 



the present calculation is based. We note |19] here that care has been taken in the above 



three-dimensional formulation to avoid double counting problems. 

In the near threhold energy region, the two nucleon energy in the final state / is very 
low and hence the NN FSI amplitude, Tj in Eq.(l), can be calculated from a number of 
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realistic A^A^ potential models in the literature. In the present work we use the A^A^ model 



developed by the Bonn group to calculate the FSI. This model is defined by a three 
dimensionally reduced BBS version of the Bethe-Salpeter equation 



T = V + ViGT 



(3) 



where G denotes the BBS two-nucleon propagator, consistent with those appearing in Eq.(0). 



(note that our definition of ViT) differs from that of Ref. |20| by a factor of —i). 

The A^A^ initial state interaction (ISI) amplitude, t/^'' in Eq.([l|), must be calculated at 
incident beam energies above 1.25 GeV . There exists no accurate A^A^ model for performing 
such a calculation. For example, the model developed in Ref. [^] can only give a very 
qualitative description of the A^A^ scattering phase shifts at energies above 1 GeV. In the 



present work we therefore follow Ref. [^| and make the on-shell approximation to evaluate 
the ISI contribution. This amounts to keeping only the 5— function part of the Green function 
Gi in evaluating the loop integration involving iG'\^'^T^^\ The required on-shell A^A^ ISI 
amplitude is obtained from Ref. [^]. As has been discussed in Ref. p2|, this is a reasonable 
approximation to the full A^A^ ISI. In this approximation, the basic effect of the A^A^ ISI is 
to reduce the magnitude of the meson production cross section. In fact, it is easy to see 
that the angle-integrated production cross section in each partial wave state i is reduced by 
a factor of [E2l 



t (r^,(p)e^2^'(^) + 1) 



(4) 



In the above equation, and //^(p) denote the phase shift and corresponding inelasticity, 
respectively; p stands for the relative momentum of the two nucleons in the initial state. 

There are a number of different approaches in the literature which model the production 
current J defined in Eq.(§) based on meson exchange models. Following Refs. []24| , p5| , we 
split the MN T-matrix of Fig.l into the pole (T^^) and non-pole (T^^) parts and calculate 
the non-pole part in the Born approximation. Then, the MN T-matrix can be written as 



■MN 



rpP I rpNP 



MN 5 



where 



-"-MN 



fldNB^dBfMNB 



(5) 



(6) 



B 



with J'mnb and qb denoting the dressed meson-nucleon-baryon {MNB) vertex and baryon 
propagator, respectively. The summation runs over the relevant baryons B. The non-pole 
part of the T-matrix is given by 



T. 



NP 
MN 



yNP , yNP-nr^P 



MN'' 



MN 5 



(7) 



where V^}^ = Vmn — Vmn^ with V^^v denoting the pole part of the full MN potential Vmn- 
Vl'fj^ is given by equation analogous to Eq.(|D with the dressed vertices and propagators 
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replaced by the corresponding bare vertices and propagators. We neglect the second term of 
Eq.(7) and hence the full MN T-matrix in Eq.(||) is approximated as Tmn — ^mat + ^mn- 
With the approximation described above, the resulting current J consists of baryonic and 
mesonic currents. The baryonic current is further divided into the nucleonic and nucleon 
resonance [N*) currents, so that the total current is written as 

J Jnuc ~)~ Jres ~l~ Jmec • (8) 

The individual currents are illustrated diagrammatically in Fig. ^. Note that they are 
all Feynman diagrams and, as such, they include both the positive- and negative-energy 
propagation of the intermediate particles. The nucleonic current is constructed consistently 
with the A^A^ potential in the BBS equation (3). The mesonic current consists of the rjpp, 
rjuju), and r^aoVr exchange contributions. The resonance current consists of the 5*11 (1535), 
Pii(1440) and Di3(1520) resonances excited via the exchange of vr, rj, p and uj mesons. 
Details of our model for the production current are given in the appendix. 



III. VECTOR MESON EXCHANGE DOMINANCE MODEL 

In order to allow for a close comparison of the model described in the previous section 
with the models |ll3HT5[] based on p meson exchange dominance, we have also constructed 



a model in which the S'ii(1535) is excited through the exchange of vector mesons. For this 
purpose, we follow Refs. to define the vector meson couplings with 5*11(1535) by using 

the following Lagrangian densites 

C^NN*{x) = -guNN*i'N*{x)'^^lt,uj^'{x)i)N{x) + h.c. , (9a) 
CpNN*{x) = -fi'pAr7v*^7v*(x)757^r ■ {3" p^ {x))^n{x) + h.c. . (9b) 

Note that the above 757^ coupling, which violates gauge invariance, is rather different from 
the tensor coupling (see Eqs.( |17c| , |17d|) ) used in our model described in the previous section 
and detailed in the appendix. 

In addition to using the 757^ coupling in the fA^ 5*11(1535) vertex (t> = p,Lj), here we 
assume the extreme case that 5*11(1535) is excited exclusively via the exchange of p and 
UJ mesons. Furthermore, for simplicity, we neglect all other resonance contributions in the 
resonance current. This is a reasonable simplification since we find that the resonance 
current contributions apart from that due to 5*11(1535) are very small. The nucleonic and 
mesonic currents are identical to the model described in the previous section. We refer to 
this model as the vector meson exchange dominance model. 

We find that this model can also roughly reproduce the total cross section data by 
choosing the coupling constants gpNN* = —0.85 and g^NN* = ~ 1.10 in Eq.d). Overall, 
the values of these coupling constants, including the signs, are consistent with those used 
in Refs. [P|,ll5|, in spite of the fact that there the ISI and FSI are treated differently from 



the present work. We therefore can use this model to investigate the differences between 
our full model described in section II and the previous works P!B|-P3| based on p exchange 
dominance. 
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IV. RESULTS 



In this section we shall present our results on the rj meson production in both the pp and 
pn collisions based on the models described in the previous sections. The parameters of the 
considered models are given explicitly in the appendix. In short, the coupling constants and 
range of form factors for meson-baryon-baryon vertices are chosen to be consistent with the 



Bonn potential in conjunction with the values used in Refs. ||2^ , |25[| and the values extracted 
from Particle Data Group [^. Thus in our calculations there is not much freedom for 
adjusting parameters. 

The total cross sections as a function of excess energy predicted by our model (described 
in section II) are shown in Fig. ^. The full results are the solid curves which are in general 
in good agreement with both the data of pp(upper panel) and p?T,(lower panel) collisions. 
For small excess energies, our pp results underestimate the data. This is usually attributed 
to the rjN FSI, which is not accounted for in the present model. Note that the results for 
pn pnrj with excess energy Q > 50 MeV , corresponding to an incident beam energy 
larger than 1.3 GeV , should be interpreted with caution, as no reliable NN phase shift 
analyses for T = states exist at present for energies above 1.3 GeV [^. To see the 
dynamical content of our model, we also show in Fig. 3 the results calculated from keeping 
only nuleonic current (dashed curves), mesonic current (dash-dotted curves), and resonance 
current (dotted curves). We see that the total cross sections are obviously dominated by the 
resonance current, and more specifically by the strong 5*11(1535) resonance (see Fig. |^). Our 
nucleonic current contributions (dashed curves) are much smaller than the resonance current 



contributions. This is rather different from the findings of Ref. |Tg]; there, instead of the 
resonance current, the shorter range part of the nucleonic current gives a large contribution 
to the pp —>■ pprj cross section. It will be interesting to know whether their model can also 
give a good description of pn —>■ pnrj data, as achieved here(lower panel of Fig.3). 

To examine the differences between our model and previous work, we show in Fig. ^ the 
results from calculations including only the 5*11(1535) resonance (solid curves) contribution. 
Within our model, this resonance excitation is due to the exchange of vr, rj, p, and u. To 
see the relative importance between these different meson exchange mechanisms, we also 
show in Fig ^ the results from vr exchange (dashed curves), rj exchange( dash-dotted curves), 
and p exchange (dotted curves). Although the u exchange is included in the calculation, its 
contribution is not shown here separately because it is much smaller than the p exchange 
contribution. As can be seen, the dominant contribution is due to vr exchange followed by 
rj exchange. The p exchange contribution is very small. Several observations are in order 
here: 



1) In contrast to the result of Refs. [|T^,0, our model, as given by Eg. ( |17cl| ), does not allow 
the 7^7^^ coupling in the pNN* vertex for the considered negative parity 5ii(1535) 
resonance. Such a coupling would prevent us from determining the pNN* coupling 
from radiative decay A^* 'jN in the vector meson dominance model (VMD) as 
explained in the appendix, since it violates the gauge invariance constraint which is an 
essential element of VMD. The simplest way of satisfying gauge invariance is to omit 
the 7^7^^ coupling and use only the tensor ^^a^^ coupling, as given in Eq.( |17d|) . This 
choice of the pNN* coupling, combined with the corresponding coupling constant as 
given in Table | - which is close to the low limit of the range determined from the 
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measured radiative decay widths - leads to a very small p exchange contribution to 
the cross section as shown in Fig. ^ (dotted curves). This is the main origin of the 
differences between our results and that of Refs. |1^,|15| where the pNN* vertex for 
A^* = S'ii(1535) is specified by the 7^7^^ coupling (see also Fig. |]). An alternative to 
avoid the gauge invariance problem while keeping the 7^7^^ term is to use a vertex 
of the form 7^(7^^^ — {itln* + mN)q^] ||29[|. Pena et al. fl^, on the other hand, have 



used the vertex 7^[7^ — {jtin* + miv)g'^/m^] in conjunction with the coupling constant 
determined from a quark model p9|. This vertex yields the same meson production 
amplitude as that of a pure 7^7^ vertex while satisfying the gauge invariance constraint, 
although only in the on-shell limit (g^ = m^). They found a significant contribution 
of p exchange to the excitation of S'ii(1535) in pp pprj. Further experimental 
information is needed to determine whether the 7^7^^ coupling is required for vector 
meson exchange. 

The 7] exchange contribution is relatively large in the present calculation. In the case 
of pp — i> ppt] its contribution to the cross section is about half of that due to the tt 
exchange. The 77 exchange contribution is subject to a relatively large uncertainty 
which arises, apart from the introduction of the phenomenological form factors, from 
the uncertainty in the rjNN coupling strength as discussed in the appendix. The 
relatively large contribution of 7] here results from using the rjNN coupling constant 
of QrjNN = 6.14, as used in the construction of the Bonn A^A^ interaction [^. This 
is close to the upper limit of the range determined empirically as mentioned in the 
appendix. However, the r] meson exchange in the Bonn potential |^ represents the 



exchange of a {J^,T) = (0^,0) quantum number and not necessarily of a genuine 
T) meson. On the other hand, the value of Qr^NN = 6.14 together with the rj — rj' 
mixing angle of 6p ~ —9.7°, as suggested by the quadratic mass formula, and the 
71 NN coupling constant of gnNN = 13.45 leads through SU(3) flavor symmetry to the 
ratio D/F ~ 1.43. This is not too far from the value of D/F = 1.73 extracted from 
a systematic analysis of semileptonic hyperon decays |3y]. Anyway, in the present 
calculation for pp ppr], the 1] exchange interferes constructively with the dominant 
TT exchange contribution, yielding the total contribution as shown by the solid line in 
Fig. ^ For pn — *■ pnr], the t] exchange interferes constructively with the tt exchange in 
the T = 1 channel (as in the case of pp — > pprj), but destructively in the T = channel 
due to the isospin factor —3 in the vr exchange amplitude. 

3) The correct description of both pp pprj and pn —>■ prirj reactions depends not only 
on the isospin factors associated with the isovector and isoscalar meson exchange but, 
also on a delicate interplay between the A^A^ FSI and ISI in each partial wave. While 
the A^A^ FSI enhances the total cross section, the A^A^ ISI has an opposite effect (see 
discussion in section II). In this connection, we mention that in Ref. |jl3| the reduction 
factor due to the A^A^ ISI is estimated to be about \(3p^) = (0.77)^ = 0.59 due to 
the state and A(ip,) = (0.73)^ = 0.53 due to ^Pi. In our calculation, however, 
the corresponding reduction factors are about X(spg) = 0.19 and A(ip^) = 0.27 near 
the threshold energy. This large discrepancy between the results of Ref. [|l3l and 
ours is due to the fact that, whereas our reduction factor is given by Eq.(^), the 
reduction factor in Ref. |jl3| is given by Aj = r]f = (e--f''"{5i))2_ argue that the latter 
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formula is inappropriate for estimating the effect of the A^A^ ISI since it exhibits a 
pathological feature: namely, when the absorption is maximum {rji = 0), this formula 
yields A, = 0, implying the total absence of the NN elastic channel and thus not 
allowing the production reaction to occur. However, scattering theory tells us that 
when the absorption cross section is maximum, the corresponding elastic cross section 
does not vanish, but is 1/4 of the absorption cross section. Note that this feature is 
present in Eq.(^. Furthermore, the authors of Ref. |TB[ apparently have identified 
incorrectly the inelasticity rji with cos^(pj), where pi is one of the two parameters (the 
other is the phase shift) given in Ref. The phase shift parametrization given in 
Ref. differs from the standard Stapp parametrization. It is obvious that with a 
more appropriate estimate of the reduction factor Aj as given by Eq.(|) the result of 
Ref. []13[ would underpredict considerably the cross section data. 



We now turn to exploring the vector meson exchange dominance model described in 
section III. This model does not have the vr and rj exchange mechanisms for exciting the 
S'ii(1353) resonance. The nucleonic and mesonic currents of our full model are kept. With 
the coupling constants QpNN* = —0.85 and Qujnn* = —1.10 in Eq.(^, we can describe both 
the pp — >■ ppt] and pn pni] data. The results are the solid curves in Fig. ^. In the same 
figures we also show the contributions from the nucleon resonance (dotted curves) and the 
nucleonic (dashed curves) and mesonic (dash-dotted curves) currents. Although the total 
cross section is underpredicted for excess energies Q > 60 MeV, it is interesting that these 



results are in line with the findings of Refs. [13-p!5| 



The results presented above indicate that the total cross section data alone cannot dis- 
tinguish two different meson exchange mechanisms for the excitation of 5*11(1535) resonance. 
It is therefore necessary to consider more exclusive observables. Fig. |^ shows the angular 
distribution(solid curves) of pp — > pprj predicted by our model (described in section II) at 
an excess energy of Q = 37 MeV. The data from Ref. ||^ are also shown. Again, the 
resonance contribution (dotted curve) dominates the cross section. As pointed out in Ref. 
T3| , the shape of the angular distribution of the latter contribution bends upwards at the 



forward and backward angles due to the tt exchange dominance in the 5*11(1535) resonance 
contribution. However, due to interference with the nucleonic (dashed) and mesonic (dash- 
dotted) currents, the shape of the resulting angular distribution (solid curve) is inverted 
with respect to that of the resonance current contribution alone. As one can see, although 
the overall magnitude is rather well reproduced, the rather strong angular dependence ex- 
hibited by the data is not reproduced by the model described in section II. At this point 
one might argue that the excitation mechanism of the 5*11(1535) resonance as given by our 
model is not correct and that, indeed, the p meson exchange is the dominant contribution. 



as has been claimed in Ref. |13|. This can be studied by considering the predictions of our 
vector meson exchange dominance model described in section HI. The angular distribution 
predicted by this model is shown in Fig. |^. Here we see that the shape of the calculated 
angular distribution(solid curve) is in better agreement with the data, although the strong 
angular dependence exhibited by the data - which shows contributions of higher partial 
waves than L = 1 - is not quite reproduced. Judging from the level of agreement between 
the two predictions and the data, one cannot discard our model in which the 5ii(1535) is 
mainly excited by vr and rj exchange in favor of the p exchange dominance model. In this 
connection, we mention that new data from COSY which will become available soon shows 
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a flat angular distribution pT 



From the above considerations, we conclude that, at present, the excitation mechanism 
of the S'ii(1535) resonance in NN collisions is still an open question. Indeed, we have just 
offered a scenario which is as good as the p exchange dominance model in reproducing the 
available data. It is therefore of special interest to seek a way to disentangle these possible 
scenarios. In this connection, spin observables may potentially help resolve this issue. As 
an example, we present in Fig. |^ the analyzing power at Q = 10 MeV (upper panel) and 
Q = 37 MeV (lower panel). The predictions of the present model are shown as the solid 
curves, whereas the predictions assuming vector meson exchange dominance for exciting the 
5*11(1535) resonance are shown as the dashed curves. The different features exhibited by the 
two scenarios for the excitation mechanism of the 511(1535) is evident. According to Ref. 
T3| , the p exchange contribution is expected to lead to an analyzing power given by 



Ay = A;^'^" sm{2e) , (10) 

where A^"'^ is positive for low excess energies, peaking at Q ~ 10 MeV and becoming 
negative for excess energies Q > 35 MeV. The corresponding results are also shown in 
Fig. 1^ as the dotted curves. Although Eq.(|TD|) gives rise to a larger analyzing power at 
Q = 10 MeV, it is interesting to see that, at Q = 37 MeV, it yields a result that nearly 
coincides with the prediction of the vector meson exchange dominance given in section III. 



V. SUMMARY 

The production of rj mesons in A^A^ collisions has been discussed within a relativistic 
meson exchange model of hadronic interactions, where the production current has been con- 
structed consistently with the A^A^ FSI used. Special emphasis has been paid to investigate 
the possible excitation mechanisms of the 5*11(1535) resonance, which is currently a subject 
of debate. It has been shown that, not only the vector meson dominance for exciting the 



5*11(1535) as advocated by Wilkin and collaborators |T^, but also the excitation mechanism 
of this resonance mainly via exchange of vr and t] mesons can describe the existing data in 
both the pp and pn collisions equally well. We have found that the analyzing power may 
offer an opportunity to disentangle these reaction mechanisms. 

A consistent description of the meson production reaction in pp and pn collisions is not 
a trivial task. As we have seen, this depends not only on the different isospin factors in the 
production current which change the relative importance of different reaction mechanisms 
from pp to pn collisions, but also on a delicate interplay between the A^A^ FSI and ISI. It 
is clear that the pp —>■ ppM and pn —>■ pnM as well as the pn —>■ dM reaction should be 
investigated in a consistent way. Also more exclusive observables than the total cross section 
such as the spin observables should be studied. 

Finally, we emphasize that the results presented in this paper should be interpreted with 
caution. The reason for this is that, as mentioned before, the A^A^ ISI is only accounted 
for using the on-shell approximation. While this may be a reasonable approximation for 
calculating cross sections, it may introduce rather large uncertainties in the calculated spin 
observables. Efforts to improve this will be published elsewhere. 
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VI. APPENDIX: PRODUCTION CURRENTS 

The T) meson production current consists of nucleonic, mesonic and resonance currents 
as shown in Eq.® and illustrated diagrammatically in Fig. ^. In the following subsections 
we construct each of these currents. A general remark to be mentioned here which applies 
to all of the currents constructed in the following subsections is that, as a consequence 
of using a three dimensional reduction of the Bethe-Salpeter equation in evaluating the 
total amplitude in Eq.(|l|), the time components of the intermediate particles involved in the 
production current suffer from an ambiguity in their definitions. In order to be consistent 
with the NN interaction used in the present work, which has been constructed by using the 
BBS three dimensional reduction, the time component of the four-momentum of a virtual 
meson at the MNN vertex is taken to be qo = e{l) — with e{l) and £:(/') denoting 

respectively the energies of the nucleon before and after the emission of the virtual meson. 
e{l) = \/P + m\j. The time component of the intermediate baryon in the nucleonic and 
resonance currents are taken to be po = ^{k) + e{p') at the B ^ M + N vertex, while at 
the N ^ M + B vertex we take Po = e{p') — uj{k). Here, uj{k) stands for the energy of the 
meson produced in the final state. 



A. The Nucleonic Current 

The nucleonic current is defined as 

June = Yl (r.-^5,f/ + ms,T,) , (11) 

with Tj denoting the rjNN vertex and Sj the nucleon (Feynman) propagator for nucleon j. 
The summation runs over the two interacting nucleons, 1 and 2. U stands for the meson- 
exchange NN potential. It is, in principle, identical to the potential V appearing in the 
NN scattering equation, except that here meson retardation effects (which are neglected in 
the potential entering in Eq.(|^)) are kept as given by the Feynman prescription. 

The structure of the t]NN vertex, Tj, in Eq.(^) is derived from the Lagrangian density 

C{x) = -g^nNN'4'N{x)-f5 (^iX + -j^l''^ i)N{.x)d^r]{x) , (12) 

where g-nNN denotes the rjNN coupling constant and A is the parameter controlling the 
pseudoscalar(ps) - pseudovector(pv) admixture. ri{x) and ipN^x) stand for the rj and nucleon 
field, respectively, and miy denotes the nucleon mass. 

The coupling constant g^,jNN is poorly known at present. The empirical values for gr^NN 



range anywhere from 1 to 7 0,^,^,0 . The values extracted from the rj photoproduction 
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analysis tend to be in the low side of this range [§3|, while a value of gr,NN = 6.14 has been 
used in the NN scattering analysis by the Bonn group [^. In the present work we use the 
value of QnNN = 6.14, consistent with the A^A^ potential V appearing in Eq.^). Also, we 
take the pure pseudovector coupling, A = 0. 

The rjNN vertex derived from Eq.(^) should be provided with an off-shell form factor. 
Following Ref . , we associate each nucleon leg with a form factor of the following form 



with An = 1.2 GeV. denotes the four-momentum squared of either the incoming or 
outgoing off-shell nucleon. We also introduce the form factor given by Eq. (p!3D at those MNN 
vertices appearing next to the r^-production vertex, where the (intermediate) nucleon and 
the exchanged mesons are off their mass shell (see Fig. Therefore, the corresponding form 
factors are given by the product F/v(p^)Fm(?^), where M stands for each of the exchanged 
mesons between the two interacting nucleons. The form factor FM{q'^) accounts for the 
exchanged meson being off-shell and is taken to be consistent with the considered Bonn A^A^ 
potential used for generating the final NN scattering wavefunction. 



B. The Resonance Current 

The production of r) mesons in A^A^ collisions is thought to occur predominantly through 
the excitation (and de-excitation) of the S'ii(1535) resonance, to which the i] meson couples 
strongly. In the present work, we also consider the Pii(1440) and Di3(1520) resonances. 

The resonance current is composed of the spin- 1/2 and spin-3/2 resonance contributions 

Jres = Jres ^ "I" '^res ^ ■ (l'^) 

The spin- 1/2 resonance current, in analogy to the nucleonic current, is written as 

Jili^^ = 5Z 5Z (^^vjN'iSN'UN* + UN'iSN'Tr^jN*^ . (15) 

j=l,2 Af* 

Here T^^jN* stands for the rjNN* vertex function involving the nucleon j. Sn*{p) = {"^ + 
niN*) / {p"^ —m%t +imN*^ N*) is the A^* resonance propagator, with niN* and Tn* denoting the 
mass and width of the resonance, respectively. The summation runs over the two interacting 
nucleons, j = I and 2, and also over the spin-1/2 resonances considered, i.e., A^* = 5*11(1535) 
and Pii(1440). In the above equation Un- (Un-) stands for the A^A^ A^A^* (A^A^* A^A^) 
meson-exchange transition potential. It is given by 

Un* = 5Z ^MNN*{q)i^M{q'^)'i^MNN{q) + ^ r^^^^.(g)zL)^^(M)(g)r^j^^(g) , (16) 

M=iT,ri M=p,u) 

where ^M{,q^) and D^y(^M){.q) denote the (Feynman) propagator of the exchanged pseu- 
doscalar and vector meson, respectively. TMNN{.q) and r^^^^(g) denote the pseudoscalar 
and vector MNN vertex, respectively. These vertices are taken consistently with the A^A^ 
potential V appearing in Eq.(^) except for the type of coupling at the ttNN vertex and 
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the uNN coupling constant. Following the discussion in Ref. |2^, we use the pv-coupling 
(A = 0) instead of the ps-coupling (A = 1) at the nNN vertex. Also, following the discussions 
in Refs. p^j2^ , |25[] , the value of the uNN coupling constant is taken to be g^jMN = 11.76. 
These exceptions apply to the currents Jres^^ and Jmec as well. Analogous expression to 
Eq.® holds for JJ^*. 

Following Ref. the transition vertices Tmnn 
resonances are obtained from the interaction Lagrangian densities 



and r^^^. in Eqs. (|T3|,|TB|) for spin-1/2 



'-riNN* 



'(±) 

-TriVAf* 



(±) 

uiNN* 



'(±) 

-pNN* 



X] 



[X 



where nix) 



TgnNN*i'N*{x) 



QujNN* 

rriN* + "^A^ 

gpNN* 

rriN* + fn]^ 
)^{x), p^{x) 



^Ar(±) + 



l-A 



rriN* ± niN 
1 - A 



rj{x) 



m^* ± rriN J ^ 
{^''uJ^'{x))'^lJN{x) + h.c. 



T ■ TC[X 



ipNix) + h.c. , 
i/^Nix) + h.c. 



fj-u \ 

x)rW(T^^f ■ {d''p^'{x))^N{x) + h.c. 



(17a) 
(17b) 
(17c) 
(17d) 



and ijjN*ix) denote the vr, cu, p and spin-1/2 nucleon resonance 
fields, respectively. The upper and lower signs refer to the even(+) and odd(-) parity reso- 
nance, respectively. The operators F^^^ and T^/f^ in the above equations are given by 



75 

r(-) = i r(-) 



757m 
7m • 



(18) 



The parameter A in Eqs.( |17a| , |17b| ) controls the admixture of the two types of couplings: ps 
(A = 1) and pv (A = 0) in the case of an even parity resonance and, scalar (A = 1) and vector 
(A = 0) in the case of an odd parity resonance. On-shell, both choices of the parameter A 
are equivalent. In this work we take A = 0. Note that we have not allowed the coupling T^j^^ 
in Eqs.( |l 7c| , |l 7d| ) in contrast to Refs. [p!4| , p!5| . Unlike the vNN vertex [v =vector meson), 
this coupling at the vNN* vertex prevents us from estimating its strength using the VMD 
because of the violation of gauge invariance. Although gauge invariant vertices which include 
the T^ij^^ coupling can be constructed |2^ , we have omitted this coupling in the present work 
for simplicity. 

Similar to the case of spin-1/2 resonances, the spin-3/2 resonance current is written as 



j(3/2) 



EE 

j=l,2 N* 



1 ^ljj^*lO, 



(19) 



Here F^^^* 



stands for the rjNN* vertex function involving the nucleon j. 



{i) mAT.) {-Qi^v + Iplvl'i + {ipVv - VvHvM'ifnj-^* + 2p^p;,/3m^.} / {y^ - + imN*TN*) 
is the spin-3/2 Rarita-Schwinger propagator. The summation runs over the two interacting 
nucleons, j = I and 2, and also over the spin-3/2 resonances considered, i.e., A^* = Di3(1520) 
in the present work. In the above equation f/^* [U'^,) stands for the A^A^ — > A^A^* {NN* — »• 
A^A^) meson-exchange transition potential. It is given by 



{q)iD\v{M){q) (q) , (20) 



M=T7,ri 



M=p,ui 
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where V^.j^n 



(g) and T^^^y^* (g) denote the pseudoscalar and vector MNN* vertex, respec- 



holds for U^* 



tively. An analogous expression to Eq.( 

The MNN* vertices involving spin-3/2 nucleon resonances in Eqs. ([T9|j20|) are obtained 



from the Lagrangian densities [32 



\ J 



(21a) 
(21b) 



2m 



^^^.(x)e^,(^)rfV7v(x)a;"^(x) 



N 



(2) 
9ujNN 



- (9A^5^*(x)e^,(^)rw^^(x)) oj^^{x) + h.c. 



(21c) 



dpNN* 

2mN 



tlj^^,{x)e^,{z)T^^^T^r^{x) ■ p'^'^ix) 



] (9,^5^.(x)e,,(z)rWf^^(a;)) ■ p'^ix) + h. 



c. 



(21d) 



where 6^,^(2;) = g^u ~ (-2 + 1/2)7^7,^. In order to reduce the number of parameters, we take 



—1/2 in the present work. 



uj^^ix] 



= d^uj''{x)-d''uj^{x) andp^'^ix) = ^ {x) - p^ ' 



X) 



The coupling constants Qmnn* used in the present work are displayed in Table.|. They 
are determined from the centroid values of the extracted decay widths (and masses) of the 
resonances from Ref. whenever available. Those involving vector mesons, are estimated 
from the corresponding radiative decay width in conjunction with the VMD. In order to 
reduce the number of free parameters, the ratio of the vNN* {v = p, u) coupling constants 
for the spin-3/2 Di3(1520) resonance has been fixed to be g^^N* / 9nNN* ~ This is 



not too far from the ratio of g^^N* / 9'^%n* ~ ~2.1 for A^* = P33(1232) extracted from 
the ratio of E2/M1 = —2.5% as determined from pion photoproduction measurements 
As for the coupling constant gvNN* {v = p, cu) for A^* = S'ii(1535), we use a value 



close to the lower limit of the range determined from the radiative decay widths given in 
Ref. ||2^ in order to emphasize the pseudoscalar meson exchange dominance in exciting 
the S'ii(1535) resonance. Since the Pii(1440) resonance is below the r]N threshold, the 
corresponding coupling constant gr^NN* has been determined by folding the results with the 
mass distribution of the resonance which is assumed to be given by a Breit-Wigner form. 
The signs of the coupling constants are chosen consistently with those used in the vr and rj 
photo-production analysis P^ , PB |. 



Following Ref. ||25| and, in complete analogy to the nucleonic current, we introduce the 
ofT-shell form factors at each vertex involved in the resonance current. We adopt the same 
form factor given by Eq.(|l3|), with mjv replaced by rriN* at the MNN* vertex, in order 
to account for the A^* resonance being off-shell. The MNN* vertex, where the exchanged 
meson is also off-shell, is multiplied by an extra form factor FAf(g^) in order to account 
for the off-shellness of this meson (see Eqs.(^^OD). The corresponding full form factor is. 
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therefore, given by the product FN{jP')FM{q^), where M stands for the exchanged meson 
between the two interacting nucleons. The form factor FM{q^) is taken consistently with 
the NN potential V in Eq.(^; the only two differences are the normalization point of 
Fv{.'f')iv = p,u!) and the cutoff parameter value of FT,{q^). Here, the form factor for vector 
mesons F^{q^) is normalized to unity at = in accordance with the kinematics at which 
the coupling constant Qvnn* was extracted, i.e., F^i^q^) = (A^/(A^ — g^))^. For the pion 
form factor F-^^q"^), following the discussion in Refs. [^,^,^], we use the cutoff value of 
= 0.9 GeV. We also use this value of the cutoff in the form factor at the ttNN vertex in 
Eqs.( [TB|pil| ) as well as in the nNN vertex appearing in the mesonic current constructed in 
the next subsection. 



C. The Mesonic Current 

For the meson-exchange current we consider the contribution from the r]vv vertex with 
V denoting either a p or u meson. This gives rise to the dominant meson-exchange current. 
The rjvv vertex required for constructing the meson-exchange current is derived from the 
Lagrangian densities 

C,U^) = -^e^p,,{d''u;^{x)Wu;^{x))7]{x) , (22) 

where Eaisufi is the Levi-Civita antisymmetric tensor with £0123 = ^1- The vector meson- 
exchange current is then given by 

{K)],zD^p{K)T%{K, k'^)zD,,{k'^)[T:^^{K)],} , (23) 

V=p,U> 

where Da/s^k^) and D^y{k'^) stand for the (Feynman) propagators of the two exchanged 
vector mesons (either the p or a; mesons as f = p or cu) with four-momentum k^ and fc^, 
respectively. The vertices involved in the above equation are self-explanatory. The vNN 
vertex V^^j^iv = p, u) is taken consistently with those in the potential used for constructing 
the A^A^ T- matrix in Eq.(|^). The uNN coupling constant is, however, taken to have the 
same value mentioned in the previous subsection. 

The coupling constant grim is determined from a systematic analysis of the radiative 
decay of pseudoscalar and vector mesons in conjunction with VMD. This is done following 
Refs. p^. 25|j37[l , with the aid of an effective Lagrangian with SU(3) flavor symmetry and 



imposition of the OZI rule The parameters of this model are the angle av{o:p), which 
measures the deviation from the vector (pseudoscalar) ideal mixing angle, and the coupling 
constant of the effective SU(3) Lagrangian. They are determined from a fit to radiative 
decay of pseudoscalar and vector mesons. The parameter values determined in this way 
in Ref. (model B), however, overpredict the measured radiative decay width of the rj' 
meson Therefore, we have readjusted slightly the value of the coupling constant of the 
SU(3) Lagrangian in order to reproduce better the measured widths. We have ay — 3.8° 
and ap = —45°, as given by the quadratic mass formula, and the coupling constant of the 
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effective SU(3) Lagrangian of G = 7 in units of Xj ^m-uvn'^^ where ra^ and stand for the 
mass of the two vector-mesons involved. The sign of the couphng constant G is consistent 
with the sign of the p'n^ and lO'n^ couphng constants taken from an analysis of the pion 
photoproduction data in the ~ 1 GeV energy region |38|. With these parameter values we 
obtain 



Qr^pp = Gcos{ap) = 4.94 



G ^-\/2 sin^(ay) sin(ap) + cos^(ay) cos(q;p) ) = 4.84 



(24) 



The r]vv vertex {y = p,uj) in Eq.(^), where the exchanged vector mesons are both off 
their mass shells, is accompanied by a form factor. Following Ref. p^ , |25[| , we assume the 
form 



F (h^ h 



/2n 



A? 



A? 



(25) 



It is normalized to unity at = ml and = 0, consistent with the kinematics at which 
the value of the coupling constant g^jw was determined. We adopt the cutoff parameter 
value of At, = 1.45 GeV as determined in Ref. |Q from the study of the u and (f) meson 
production in pp collisions. This form factor has been also used in the study of the rj' meson 
production in Ref. p5| . 

Another potential candidate for mesonic current is the ?7ao7r-exchange current, whose 
coupling constant may be estimated from the decay width of Oq into an t] and vr. We take 
the Lagrangian density 



^rjaon ('^) 



(26) 



where ao{x) stand for the Oo meson field. m,r and rua^ stand for the masses of the vr and 
tto meson, respectively. Using the measured decay width from [^] we obtain a value of 
~ 1.81. The sign of this coupling constant is not fixed. We assume it to be positive in 
the present work. Since the contribution of the rjaoTc current is small, the sign of its coupling 
constant will not affect the major conclusion of the present work. 
The 7]aoTT current reads 

Jrja^n = { [F^^^^] (fc^ JR^^^^ ( A;^^ , k^)iA^{kl)[T N {k^)h} + (1^2) (27) 

in our previously-defined notation. The vertex Fa^jVAr is taken consistently with that in the 
NN potential V in Eq.(^, while for the ttNN vertex we use the same one (including the cut- 
off parameter value) mentioned in the previous subsection. The t](1qTI vertex, F^fjo7r(^ao) k^^, 
is assumed to have a form factor given by 



F (P k 



do 



kl 



mz 



Al - k^ 



(28) 



with A„„ = A^ = 1.45 GeV. 
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The total mesonic current is then given by 



Jmec ^ ^ Jt]VV ~I~ Ji 



(29) 



V=p,Ul 



There are, of course, other possible mesonic currents, such as the rjcjcj)- and r/^^-exchange 
currents, that contribute to rj meson production in NN collisions. Their contributions have 



been estimated in a systematic way following Ref. [24| and were found to be negligible 
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TABLES 

TABLE I. Coupling constants qmnn* used in the resonance current. Those in parenthesis for 
vector mesons refer to {g^MNN*-'9^MNN*)- '^^^ ps-pv mixing parameter is fixed to be A = and 
the off-shell parameter z = —1/2. The masses and widths of the resonances, rriN* and Fjv*, are in 
units of MeV . 



N* 


^ii(1535) 


Pn(1440) 


Di3(1520) 




(1535,150) 


(1440,350) 


(1520,120) 


71 


1.25 


6.54 


1.55 


V 


2.02 


0.49 


6.30 


p 


-0.65 


-0.57 


(6.0, -2.1) 




-0.72 


-0.37 


(-2.1, 0.7) 
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FIG. 1. Amplitude for the NN NNr] reaction considered in the present work. Tmn denotes 
the MN T-matrix. ISI and FSI stand for the initial and final state NN interaction, respectively. 
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FIG. 2. T] meson production currents included in the present study. Upper row: nu- 
cleonic current Jnuc, M = 7r,r], p,uj,a,ao- Middle row: nucleon resonance current Jres, 
N* = 5ii(1535), Pii(1440) and Di3(1520), M = Tr,r),p,uj. Lower row: mesonic current Jmes, 

V = p,LO. 
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FIG. 3. Total cross sections for the pp — > ppr] (upper panel) and pn — > pnr] (lower panel) 
reactions as a function of excess energy within our model described in section II. The dashed 
curves correspond to the nucleonic current contribution while the dash-dotted curves to the mesonic 
current contribution; the dotted curves represent the resonance current contribution. The solid 
curves are the total contribution. The data are from Refs. ill-oJ 
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FIG. 4. Same as Fig. ^, except that it shows the Sii(1535) resonance contribution only. The 
dashed curves correspond to the vr exchange contribution while the dash-dotted curves to the r] 
exchange contribution; the dotted curves represent the p exchange contribution. The solid curves 
show the total contribution. 
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FIG. 5. Same as Fig. I but using the vector meson exchange dominance described in section III. 
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FIG. 6. Angular distribution of the emitted r] meson in the cm. frame of the total system at an 
excess energy of Q = 37 MeV. The dashed curve corresponds to the nucleonic current contribution 
while the dash-dotted curve to the mesonic current contribution; the dotted curves represent the 
resonance current contribution. The solid curve show the total contribution. The data are from 
Ref. §. 
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FIG. 7. Same as Fig. I but using the vector meson exchange dominance of section III. 
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FIG. 8. Analyzing power for the reaction pp — > ppr] as a function of emission-angle of rj in the 
cm. frame of the total system at an excess energy of Q = 10 MeV (upper panel) and Q = 37 MeV 
(lower panel). The dotted curve corresponds to the case of p exchange dominance according to 
Ref. The solid curve corresponds to the full model calculation described in section II. The 

dashed curve is the prediction assuming vector meson exchange dominance as described in section 
III. 
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